A Criterion for the Critical Number of Fermions and Chiral Symmetry Breaking in 

Anisotropic QED 3 
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By analyzing the strength of a photon-fermion coupling using basic scattering processes we calcu- 
late the effect of a velocity anisotropy on the critical number of fermions at which mass is dynamically 
generated in planar QED. This gives a quantitative criterion which can be used to locate a quantum 
critical point at which fermions are gapped and confined out of the physical spectrum in a phase di- 
agram of various condensed matter systems. We also discuss the mechanism of relativity restoration 
within the symmetric, quantum-critical phase of the theory. 
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I. INTRODUCTION 

Quantum Electrodynamics (QED), the quantum the- 
ory of radiation and its interaction with matter, is a sub- 
ject with far reaching impact in physics, from the cal- 
culation of the magnetic moment of the electron 1 - to the 
widespread use of diagrammatic tools in every corner of 
theoretical physics. However, despite this long tradition 
as a central paradigm in physics, QED continues to con- 
front us with many new challenges as its modern rein- 
carnations emerge as effective theories of strongly corre- 
lated many-body problems, from quantum spin liquids to 
high temperature superconductivity to graphene. Back 
in Feynman's day, we would have been very surprised if 
someone were trying to solve a relativistic problem with 
two speeds of light. These days, such problems are actu- 
ally ubiquitous in modern theoretical physics, one exam- 
ple being the effective theory of low energy excitations 
in a correlated d-wave superconductor. That theory is 
equivalent to an anisotropic quantum electrodynamics in 
2 + 1 dimensions 2,3 , in which different "speeds of light" 
appear naturally. Furthermore, there are other systems 
in which the low energy description reduces to different 
versions QED3, such as various forms of quantum spin 
liquids 4 -^ or the physics of graphene layers^. All of these 
different problems share a common feature: they have 
a nodal structure that resembles a relativistic spectrum 
for the low energy excitations. However, a distinctive 
feature of the QED3 particularly relevant for cuprate 
superconductors is that it contains a significant intrinsic 
anisotropy, which exists due to the difference between the 
Fermi and gap velocities (vp, va)- 

It is known that the fermionic anisotropy, a — vf/va, 
is irrelevant in the perturbative renormalization group 
(RG) sensed, as long as the system is in the symmetric 
phase of QED3. This is the quantum-critical phase of the 
theory, in which strongly interacting massless fermions 
acquire anomalous power-law behaviors in their various 
correlation functions. The anomalous dimension expo- 
nents characterizing this unusual state are universal and 
typically depend only on the total number of fermion fla- 
vors, N. This however is true only as long as N > N c , 
where N c is the critical number of fermions at which the 



fermion mass is dynamically generated via the mecha- 
nism of spontaneous chiral symmetry breaking (CSB). 
Once CSB takes place, the fermions are gapped and con- 
fined out of the physical spectrum. This heralds a differ- 
ent, massive phase of the theory which typically trans- 
lates to a different state in the underlying condensed mat- 
ter system. The results of Ref. 7 are generally valid for 
arbitrary anisotropy as long as the number of fermions 
N is sufficiently large or for small anisotropy when N is 
greater than N c of the isotropic case. 

Evidently, N c is an important number within the the- 
ory, not in the least because antiferromagnetic order in 
effective theories of high temperature superconductors 
and quantum spin liquids generically arises via the above 
non-perturbative phenomenon of CSB - for example, in 
the context of cuprates, the chiral mass generation cor- 
responds to the onset of a spin density-wave order from 
within a quantum disordered d-wave superconductor^ 
while it describes the formation of the Neel antiferromag- 
netic state and a whole family of other order states in the 
context of quantum spin liquids!^ Unlike the exponents 
of the critical massless phase, however, N c itself is not 
universal. Consequently, an important question arises 
within the anisotropic QED3: to what extent is the crit- 
ical number of fermions, N c , which defines the boundary 
between broken and unbroken chiral symmetry, affected 
when such anisotropy is present? 

In this paper, our goal is to provide an answer to this 
question. Of course, ours is not the answer, for two rea- 
sons. First, since QED3 is a strongly interacting theory 
its exact behavior is beyond our reach. Second, since N c 
is not universal there are actually many different 2V c 's: 
the nominally irrelevant couplings within the QED^ the- 
ory of a quantum disordered d-wave superconductor 8,9 
are very different than those of lattice-based quantum 
spin liquids!^ Furthermore, both these N c 'b are dif- 
ferent from the intrinsic N c of the pure QED3 field 
theory considered here (defined through Balaban-Jaffe 
regularizationj 1 * for example). However, all these issues 
notwithstanding, even in the absence of the exact solu- 
tion, it is still possible to make a rather accurate de- 
termination of the parametric dependence of N c on the 
anisotropy, once an "exact" N c is known for the isotropic 
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case from a different source, say from numerical simula- 
tions. Our goal is to devise a criterion for determination 
of N c within the anisotropic QED3 which, while not ex- 
act, still provides a rather accurate description of how 
the CSB boundary changes as a function of the parame- 
ters of the theory. The philosophy here is similar to the 
one behind the ubiquitous Lindemann^ criterion, orig- 
inally proposed to predict the melting point of a solid. 
While not exact, the Lindemann criterion has proven it- 
self a remarkably accurate and useful in a wide range of 
situations, from classical to quantum solids, from Wigner 
crystals to Abrikosov vortex lattices. 

To devise our criterion, we point out that mass gen- 
eration is a consequence of the fermion-photon coupling 
strength in QED. This strength is generally renormal- 
ized from its bare value by virtual polarizability of the 
vacuum. Relying on this fact, we propose a natural way 
to measure the strength of the gauge field by focusing on 
the matrix element that represents processes in which one 
photon is exchanged between two fermions. We stipulate 
that the CSB and mass generation take place when this 
matrix clement exceeds certain critical value. Within the 
isotropic QED3 this is manifestly an exact statement - 
the only unknown is the actual value of N c which we can 
either infer from a separate calculation or borrow from 
numerical simulations . 12 i 13 i 14 Once our "Lindemann cri- 
terion" is calibrated in this fashion, we proceed to eval- 
uate the appropriate matrix element in the anisotropic 
theory and propose that the CSB takes place when this 
matrix element exceeds the same critical value. 

Following the above procedure, our criterion allows us 
to derive the following main results: first, we show that 
N c is not just a function of the anisotropy a, but also a 
non-monotonic function of vp and va themselves. The 
reason behind this peculiar behavior is the existence of a 
third relevant velocity in the theory, the speed of light, 
c s , which naturally appears in the Maxwellian action for 
the gauge field. We confirm our findings for va > c s by 
using the Schwinger-Dyson (SD) equation within the Pis- 
arski approximation; 15 both results show the same func- 
tional dependence on va in this regime. Another result 
is that the critical number of fermions decreases even for 
a = 1, as long as vp — va^> c s . Thus, in the context of 
the QED 3 effective theory of underdoped cuprates, we 
generically expect a non-superconducting, non-magnetic 
pseudogap state to intrude prior to the emergence of anti- 
ferromagnetism as d-wave superconductor is underdoped 
toward half-filling. 

II. CRITERION AND CALCULATION 

We start by defining the form of the anisotropic QED3 
Lagrangian and the main physical quantities that will be 
used in our calculations. The Lagrangian is 

JV 

CAniso = y^^iv^lnjidn ~ a^ipi + C[a^(x)] (1) 
1=1 



where ipi are four component spinors and is a gauge 
field. In condensed matter environment, Eq. ([1]) is 
the low energy effective theory of some strongly cor- 
related electron system. For example, in a phase- 
disordered d-wave superconductor, the gauge field a M 
couples nodal Bogoliubov-deGennes-Dirac fermions to 
fluctuating quantum vortex-antivortex pairs through a 
Maxwellian action C[a^(x)\. 7^ are the standard gamma 

matrices, such that {7^,7^} = Sfiv, and vjp is defined as 
(1, Vp, va) and (1, va, vf) for nodes 2 = 1,2 respectively. 
Fermi and "gap" velocities, vp and va, define fermionic 
anisotropy at each node. 

To make further progress we first focus our attention 
on a particularly simple case: the non-interacting one 
(a M = 0). This case is equivalent to saying that the 
phase is rigid or that the departure point in our analysis 
is the fully-ordered superconducting state. In that case 
Eq. ((T|) is reduced to 

2 

Cmrac = + m )^l (2) 

1=1 

where a mass term has been introduced only for normal- 
ization purposes. In the end our results will be robust in 
the m — > limit. 

The equations of motion for these non-interacting 
Dirac fermions are: 

(-v#V^m + m )^i = ( 3 ) 

with I = 1,2. From this expression it is clear that the 
task of solving Eq. (3) can be performed independently 
for each fermion flavor 1 and 2. Now, if we focus on I = 1, 
we can reduce the anisotropic equation of motion to the 
isotropic one using a simple change of variables r' = r, 
x' = vpx, and y' = vaTJ- It is worth remarking that the 
intrinsic fermionic anisotropy cannot be removed in the 
case of the real d-wave superconductor once the gauge 
field is present, as can be seen from the definition of V/p. 

We now outline how a typical perturbative approach 
for computing the effect of coupling to the gauge field 
works: starting from the solutions obtained for the 
problem in which a M = 0, namely ipi(x) and ip2{x) 
- both having the anisotropy encoded in their spa- 
tial oscillations - we will perturbatively introduce the 
field a M . First, we choose our free solutions for each 
node in the form ipf(x) = s r (m, 0) exp(— imt) and 
ipf(x) = t r (m, 0) exp(-Mmi), where for convenience we 
have defined s x (m,0) = (1,0,0,0), s 2 (m,0) = (0,1,0,0), 
h(m,0) = (0,0,1,0), i 2 (m,0) = (0,0,0,1). The eigen- 
states of Eq. (3) are then 

(4) 

where b r p and d r $ are standard fermionic annihilation 
operators, while p and r are the momentum and spinor in- 
dices, respectively. In the case of decoupled nodes, when 
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dp = 0, all the information about the anisotropy is en- 
coded in the exponential factor. Now, for finite momen- 
tum we can use the ansatz tp^(x) — s r (fc) exp(— ik ■ x') 
and il>i{x) = t r {k) exp(+ifc-x'), where x! — (r, vfx, vav) 
for node 1. The spinors s r (k) and t r {k) can be found di- 
rectly from the equation of motion 

(#-m)s r (fe)=0 + m)t r (k) = (5) 

where If, — 7 Al fcp. Thus, applying Qt + to) and (If — m) 
to s r (m, 0) and t r (m, 0) respectively we can generate the 
needed solutions. 

Using these free fermions in the Heisenberg picture, 
we now introduce the gauge field by employing the well- 
known formal solution of the scattering matrix 5, equiv- 
alent to: 

S = T[exp(-i J d 2+1 xHi{x))] (6) 

where T and : are the time ordered operator and the 
normal order operation, respectively. The interaction 
Hamiltonian density that couples fermion fields to the 
gauge field is given by 

2 

Hi{x)=Y,-^h^i: (7) 
i=i 

Now, as usual, we want to analyze the matrix element 
< f\S\i >, where both the final and initial state are 
assumed to be eigenstates of the unperturbed Hamilto- 
nian. Using the Heisenberg representation of the free 

(a) > (b) 




FIG. 1: (a) Diagramatic form of the Schwinger-Dyson equa- 
tion at one loop approximation, (b) Diagram used to measure 
the strength of a photon. 

fields ip(x) and the Wick's theorem, we recognize that 
the only relevant diagrams are the interactions photon- 
fermion and fermion-fermion, the so called direct scatter- 
ing and the Moller scattering. The calculation simplifies 
here because, due to its fluctuating nature, the contri- 
bution of a M to the first term of the S matrix expansion 
forces the said matrix element to vanish. Consequently, 
the leading contribution relevant for our purposes comes 
from the second order term where the contraction of 
a^ixi) with a v (x2) is nothing else but the photon prop- 
agator D^^xi — X2), which is the real space inverse of 
the polarization function 

n -V N Jka {n) kJa M ^pk P 9 ( ^k x \ 



where the diagonal nodal "metric" gjZ is defined^ by 

(1) (2) 1 (1) (2) 2 i (1) (2) 2 

fltio = #00 = ! > £U = .922 = V F and 522 = 011 = V A ' 

and 7^ = \J gffl 7^ ■ We will now demonstrate that 
this second order element of the perturbation theory 
can be used to gainfully define the "strength" of the 
fcrmion-photon interaction and determine the critical 
number of fermionic flavors N c (a), without ever solving 
the Schwinger-Dyson equation, even though a perturba- 
tive approach itself is condemned to failure. 

In order to accomplish this, let us focus on the mass 
generation problem. It is well known that Schwinger- 
Dyson equation, Fig. 1(a), has a non trivial solution 
(rn > 0) once the initially soft photon (factor 1 /N in the 
large N limit) becomes harder at a critical number of 
fermions, N c ; with 3 < N c < 5^^^^. That means for 
N < N c mass will be dynamically generated and there- 
fore chiral symmetry would be broken. 

The existence of anisotropy makes to analyzing the 
birth of a gapped state through the mechanism of 
mass generation using the Schwinger-Dyson formalism 
almost impossible, with the noteworthy exception of the 
isotropic limit and the so called small anisotropy case 
(vf = 1 + S,va = 1, where S <C l)l i 18 i 19 . Our aim is to 
explore mass generation for arbitrary anisotropy, because 
this is the physically relevant regime. Our technique is 
based on the observation that as we lower N we are ef- 
fectively making the photon interaction stronger. This is 
clear because the photon propagator contains a screen- 
ing factor jj,. From this simple observation we can see 
that mass generation is a phenomenon that is intrinsi- 
cally tied to the strength of the gauge field (photons). 
Put simply, we can translate the meaning of N c into the 




FIG. 2: The calculated ratio f3. The minimum value of v^and 
va is 0.15 and the maximum 1.5. This graph shows that 7V C 
is not a simple function of a but a more complex function of 
Vf and «a ■ 
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following statement, which clearly is a necessary condi- 
tion for mass generation: If the photon- exchange inter- 
action is stronger than some threshold value A c then the 
fermion mass will be dynamically generated. This asser- 
tion clearly lacks the precision of a mathematical theo- 
rem but we now proceed to demonstrate that it has the 
practical virtue of a useful criterion. 

The problem that we are facing is straightforward: 
How can we usefully quantify the strength of the gauge 
field? The simplest way to measure the strength of 
a photon is to analyze a process in which one photon 
is exchanged. For the sake of simplicity, we will ana- 
lyze the fermion — fermion scattering mediated by one 
photon, as shown in Fig. 1(b) which corresponds to 
the second order term in perturbation theory. Thus us- 
ing the N c borrowed from the isotropic limit, we define 
A c =< i\S^\f > [N c (a = 1)], the matrix element of the 
S matrix evaluated at the isotropic limit, as the natural 
candidate to measure the strength of the photon. We 
can safely say that in an anisotropic theory mass will be 
generated at the critical strength defined by A c . 

Thus solving the equation, A c =< i\S^\f > [N c (a)] 
for N c (a), we will find the influence of anisotropy on 
the critical number of fermion flavors. In the isotropic 
case, the choice of node will not make any difference, 
however in the anisotropic case the situation will be quite 
different. Regarding the symmetries of the problem we 
must analyze the invariant, < i\S^\f > (vf,v a )+ < 
i\S^\f > (v A ,vp) because nodal particles can be born 
indifferently in nodes 1 or 2. This quantity is invariant 
under the transformation vp <-> v A , which is a physical 



symmetry of this problem. 

We now proceed to compute the first non zero matrix 
element < i\S^\f >. The second order process following 
from the S-matrix term, Eq. (6) is: 



5 (2) oc T[J Hi(2i)Hi{2i)dxnix2] 



(9) 



which can be written (integration with respect to x\ and 
x 2 is implicit) as a sum of terms of the type 



T[: V'lOcim V (£1)^1(2:1) :: ipi(x 2 )'yi, ' a l/ (x 2 )i)i(x 2 ) :] 

, (10) 

By applying Wick's theorem to the previous expression 
terms with and without contractions are obtained. How- 
ever, due to the fluctuating nature of the gauge field most 
of those terms will vanish. The survivor will be the one 
containing the contraction of the gauge field with itself. 
Thus, the problem is reduced to compute the matrix el- 
ement: 



: ^1(^1)^1(^2)7^ 'i>i{x2)D llv {x2 - ari) : 

(11) 

To evaluate such matrix element we must define the ini- 
tial and final states of the system \i >,\f >. As the 
strength of the photon should be independent of the 
states that we use to measure it, we are free to use 

It >= bl - bl - |0 > and I f >= b\ - b\ - |0 > (12) 

I npi r 2 p2 1 '•' r' 1 p' 1 r' 2 p' 2 ' v ' 

which we already know from the free theory The relevant 
matrix element is 



J 



1 < i\SW\f > (A,v F ,v A ) ~ — f d 3 k (—J——) Q - KpAD^g^^Q^puPi 
N v A v F J \Pio(Pio ~ k )J 



k) 



(13) 



r 



where Q^Pjff) = s(P)^' s(p), with s — s"/ - The in- 
tegral is to be performed over all the allowed k — space, 
but we have one constraint that is the spectrum of the 

Dirac fermions, k = \J (vFk x ) 2 + {v A ky) 2 , which allows 
us to further reduce the evaluation of this integral. 

The integration limits are fixed from the condition that 
the momentum transferred in the processes shown in Fig. 
1 (b) can not be larger that the momentum of the incident 



particle, pi. We have set p x = p y = A as an upper cut 
off. We have also verified that the only effect of choosing 
a different cut off for different spatial directions is that it 
makes the convergence slower. Thus, using the invariant 
measure of the strength of the gauge field, it is straight- 
forward to show that the ratio between critical number of 
fermions at velocities v F and v A respect to the isotropic 
point is 



(3 



< i\S^\f > (A,v F ,v A )+ < i\S^\f > (Wf) 
2<i|S( 2 )|/>(A,l,l) 

I 



(14) 



Thus, the critical number of fermions for the anisotropic theory is given by the product of, /3(vf, v a ), and the crit- 
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ical number of fermions at the isotropic point N c (1, 1). 
In principle, this measure of the strength of the gauge 
field depends on the initial momentum of the particles, 
but it turns out that the ratio between the strength of 
the gauge field with anisotropy and the strength of the 
gauge field without anisotropy a = 1 is not sensitive to 
the momentum of the incident fermions, pi. In fact, all 
the calculated quantities converge to a fixed value for 



a fairly low momentum cut off, A. It is important to 
remark that (3 is gauge invariant due to the transverse 
nature of the photon. 

Indeed we can do better. By isolating the leading di- 
vergent behavior of the scattering amplitude we have ob- 
tained that in the A — > oo at the leading order in the 
large N approximation, the ratio /3 is given by: 



P(v f ,va) = Ci 



'l + v 2 F E 



'A 



+ C2 



K 



K 



J A 



(15) 



where 

E(m) = [ 2 [1 - msm(8)} 1/2 d8 
Jo 

K(m) = [ 2 [l-?7isin(60r 1/2 d6l (16) 
Jo 

= ( l-2(vl + vi)-v* F vl+v% + vi \ 

= 2 - 3 (4 + vj) + (v 2 F -vlf+ v 4 f v 2 a + yjvF 
^ 3V2irv F vi^(l + v F )(l + vl) 

are the elliptic functions of first and second kind respec- 
tively, and the prefactors are polynomials that are sym- 
metric under vf *-> va- Thus (3 is explicitly invariant 
under the transformation vp <-> va, as it must be. 

An important case is the one in which there is no 
fermionic anisotropy, i.e. v F = va- In this case all 
the integrals can be performed and the full divergent be- 
havior of Eq. (p~3|) can be isolated by going to cylindrical 
coordinates as the x <-> y symmetry suggests. The re- 
sult is just the limit behavior of Eq. (|15p when taking 
vf = va- In that limit the dependence on the fermionic 
velocities becomes particularly simple. 

"-■fri (17) 

To test the accuracy of our technique, the result is com- 
pared with the exact numerical integration in Fig. 3, 
showing excellent agreement. For low speeds va < 1, 
our result suggest that the critical number of fermions 
increases. We discuss this case in more detail below. On 
the opposite limit when the gauge field velocity is small 
compared with the other two velocities the critical num- 
ber of fermions goes to zero. 

A non trivial dependence on the gauge field velocity 
was also found; there are three different regimes depend- 
ing upon the value of vp, with respect to c s , as shown in 
Fig. 2. For v F < c s , the critical number of fermions in- 
creases, reaching unexpected values for large anisotropy, 




1 2 3 4 5 



V F 

FIG. 3: Comparison between the numerical and analytical 
calculation of /3 in the case for which vf — va- 



leading to a gapped state in this sector. On the other 
hand for vf > c s , N c decreases and we can generically 
expect gapless excitations. The most unexpected result 
occurs if vf = c s , in this case N c (a) w N c (l) for a > 1. 
Even though we found small deviations with respect to 
the isotropic value near to a = 2, they do not change the 
integer part of N c (See also Fig. 5). Thus, in this regime 
large anisotropy will be completely irrelevant. 

Using AT C (1) = 4, which is the integer closest to 
the gauge invariant critical number of fermions 128/37T 2 
found by Nashi£, we have obtained a phase diagram in 
the (vf, va) space, Fig. 4. In this plot it is clear that 
if both velocities are smaller than c s , the critical number 
of fermions never goes to zero. However, outside that 
square region N c — > for high anisotropics. 

To show how important it is to consider the scale de- 
fined by c s we have also plotted /3(a) for different values 
of Vf (0.5,1,5.0) that we kept fixed as va was varied 
(See Fig. 5). In this plot we make evident that c s is a 
relevant parameter which we are allowed to set to 1, how- 
ever the detailed behavior of the system is not a simple 
function of the anisotropy, a, but an explicit function of 
both parameters vf and va- 

In order to test our results we compare them with 
recent numerical work that has been done by Thomas 
and Hands (T-H)i2 in an attempt to understand how 
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0.025 0.25 0.5 0.75 1 1.25 1 



FIG. 4: Phase diagram for the critical number of fermions in 
the velocity space. We have set the critical number of fermions 
for the isotropic theory at N c = 4. Inside the square defined 
by the lines Vf = 1 and ka = 1 the critical number of fermions 
never vanishes. The blue line corresponds to Vf = 1/va and 
represents the one dimensional domain simulated by Thomas 
and Hands^. 



anisotropy can modify the properties of QED3. The 
lattice (Euclidean) version used by T-H includes an 
anisotropy which is intended to mimic the contin- 
uum model in the incarnation presented by Lee and 
Herbuti^. The QED% theory presented in Lee's article is 
parametrized in terms of two quantities 8 — y/vpVA and 
A = vf/va- T-H have used an extended lattice model 
similar to the one used by Dagotto et al2£, that in the 
continuum limit resembles the behavior of Eq. (1). In 
order to perform the simulation the action used was: 



5 = E E fl 8 Xi(sR,,'»M + ~a 8 e£„0»0 (is) 

i— 1 x : x' 

where the anisotropy was introduced in the fermion ma- 
trix 



m x . x , = ±j2u x ) 



2a 



and 



(_!)*!+. 



(19) 

(20) 
(21) 



with x\ = x, X2 — y and X3 — t, is the Kawomoto-Smit 
phase of the staggered fermion field. The lattice spacing 
is a. An important definition is that of the anisotropy 



factors, A, 



A, 



and A T 



This def- 



that T-H lattice theory does not keep the flavor symme- 
try of the model relevant for cuprates^ i 3 ' 7 ' 12 . Regardless 
this intrinsic drawback of the method we will show that 
it still mimic the qualitative behavior of the flavor sym- 
metric anisotropic QED3 at least in the low anisotropy 
limit. 

In their numerical simulation T-H have a single pa- 
rameter, which is equivalent to a — vf/va, and as an 
extra constraint they have set 8 = 1. This choice imply 
that vf — 1/ua an d therefore a — l/v A . That means 
that T-H have simulated a 1-D domain of the whole pa- 
rameter space (vf,va)- That domain is shown in Fig. 
4 as a blue line. In order to make a link between the 
amount of condensate < ijjtp > and its relation with 
N c (a) we will assume that the functional form of the 
dynamically generated mass does not chang e - 16 i as we 
introduce anisotropy in the system. This can be explic- 
itly checked in the small anisotropy limits. As long as 
we are inside to the broken phase and near to the bound- 
ary between broken and unbroken phases, the mass has 
the following functional form 



-2tt 



m = m exp 



N c (v F ,v & ) 
N 



(22) 



1 



inition is important for our purposes because it shows 



where m may depend on N. However, what is impor- 
tant is how to define the boundary between the massless 
phase and the massive phase. This is done by solving Eq. 
(122p when m = 0. We can also interpret this equation 
in the following way: For a fixed number of Fermions, 
A, this equation allow us to understand how the mass 
changes as a function of N c (a) when A is close to iV c . 
In fact, for any number of fermions when N c (a) = N it 
is prohibited to have any condensate. That means that 
the critical anisotropy a c that solve m(a c ) = for any 
value of N is the same that the critical anisotropy that 
solves N c (a c ) = N. Hands reported that such a critical 
anisotropy a c 4 in a 16 3 sites lattice simulation which 
is in agreement with the critical value found by us a c ~ 3. 
This decreasing behavior is in contrast to the one found 
in Ref. 18, where it was claimed that A c increases as 
a function of the bare anisotropy. Taking into account 
that the lattice simulation is a non perturbative method 
that does not relay in any educated ansatz, T-H results 
strongly support our view of the phenomenon. Still, we 
should mention that the critical anisotropy calculated by 
T-H using the anisotropic scaling is not, in the strict 
sense, quantitatively accurate in the context of cuprates, 
for two reasons. First, the scaling used brakes the crys- 
tal isotropy, or, in different words, their simulation is not 
invariant under flavor exchange. Second, they assumed 
that the two fermionic velocities change little around the 
gauge field velocity c s = 1. Due to the collective nature 
of phase defects we expect that c s < f)p in cuprates^. 

On the other hand lattice simulations are always per- 
formed in finite lattices and therefore the correct com- 
parison with our results should be done by considering 
both, an upper cut off and a lower cut off. In principle we 
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are free to make the upper cutoff as large as we want but 
the lower cut off dependence may be important when 
comparing our results with finite lattice simulations 22 . 
Setting the upper cut off as A u — 1 then the lower cut 
off should be ~ 1/Z, where L is the size of the sys- 
tem. As we change the lower cut off we have found no 
significant differences in our results as it goes to zero. 
We have also investigated the effect of the breaking of 



1.5 













I 1 


i , 


, i , 


v F -0.5 

v F -1.0 

v r -2.0 




0.5 


1 


1.5 


2 2.5 





a 



ca 




a 

FIG. 6: /3 within the domain defined by the blue line of Fig. 
4. The continuous line represents the flavor symmetric cal- 
culation. The dotted line is the non-symmetric calculation. 
For the latter case we have only plotted the real part because 
the breaking of the flavor symmetry introduces a small but 
persistent imaginary part. 



FIG. 5: Ratio between the critical number of fermions in 
the anisotropic theory and its equivalent within the isotropic 
theory as a function of anisotropy, /3(a). The three curves 
represent the three regimes found: N c (a) < N c (l) if vf > 1, 
N c (a)/N c (l) approaches a constant ~ 1 if vf = 1, and N c (a) 
increases for increasing a if vf < 1 • 

the flavor symmetry in our scheme. For this purpose we 
have used only one amplitude for the photon, so that the 
final kernel is not invariant under vf ^a- I n this case, 
we have obtained the same qualitative behavior as shown 
in Fig. 6, and the small negative slope at the isotropic 
point, detail that is in agreement with the results shown 
by T-H in [Fig. 5, (Ref. 11)]. This clearly is a symptom 
of the breaking of flavor symmetry. However, there is 
a very important issue that emerges once we arbitrarily 
break flavor symmetry. The function f3 starts to pick up 
a phase which is unphysical. (3 must be a real number as 
N c is. This makes evident that even with the simplest in- 
teraction between fermions the flavor symmetry is needed 
in order to obtain a meaningful value of (3. This obser- 
vation suggest that numerical simulations that preserve 
flavor symmetry are the only reliable way to extract ac- 
curate critical values. However, we must encore that the 
qualitative behavior of T-H simulation agrees with the 
physical picture proposed in this article. Looking at the 
obtained phase diagram, Fig. 4, the question that nat- 
urally arises is, which regime is the physically relevant 
one? Clearly our answer will depend on the ratio vf/c s . 
At finite but still small temperature, T, we can use the 
continuum vortex-antivortex Coulomb plasma model in 
order to estimate the gauge field velocity. Identifying the 
speed of light from the Maxwellian form of the action for 
the gauge field we find that c s ~ yJni/T at finite T, 
where ni is the density of vortices^. Thus for any T > 
as we approach the superconducting state n; — > and 
thus c s — ► 0, resulting in a protected symmetric phase. 
As ni increases, for small but finite T, c s may reach very 



high values, but if those values are larger than Vf is un- 
clear. If we go to the line T = quantum fluctuations 
will drive the system into a region in which the value 
of c s will depend on the specific value of the dynamical 
critical exponent z, that in some simplified calculations 
was adopted to be 1. The more striking problem about 
identifying the precise value of c s is that this velocity is a 
function of the correlation length, c s = f(£, sc ), but at the 
same time we know that £ sc = £, S c(x, T), thus a self con- 
sistent treatment or knowledge of the correlation length 
from experiments will be needed to settle this problem 
and give an accurate phase diagram that identify the rel- 
evance of these different regimes. 

To provide further evidence of our findings we have 
also re-analyzed from a different perspective the case in 
which fermionic anisotropy is one. We have applied Pis- 
arski's technique^ to find the qualitative behavior of the 
dynamically generated mass as a function of va- 

We will assume that E (p) ~ S (0). Thus in the ap- 
propriate integration interval (E (0) , A) the mass will be 
a constant. This assumption is certainly incorrect, as 
was shown by Appelquis t 16 i 17 et al. However, it will al- 
low us to compare the qualitative behavior of the mass 
as a function of va, for the case in which vf = va- In 
this case the sums can be performed analytically by go- 
ing to cylindrical coordinates instead of spherical ones. 
The integration was performed over the shell defined by 
a lower cutoff E (0) and an upper cut off A. Thus the 
Schwinger-Dyson equation: 



d 3 k (23) 



can be solved in this rough approximation and the result 
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is: 

which is a real number as long as va > 1- Thus as long as 
the neutrinos move faster than the photons the generated 
mass indeed does depend on va- 

We expect that - given the fact the decay factor in 
Pisarski's result is of order N c - the correction found for 
that factor will give us the functional behavior of N c (va)- 
Thus, as va ^> 1 we would expect that N c ~ 1/i>a< This 
is consistent with the result obtained from the proposed 
criterion, Eq. (|17|) . However, we should warn the reader 
that even thought a mathematical expression can be ob- 
tained for va < 1 the nature of the system will change 
in this case, casting doubt on the validity of our crite- 
rion in that region. Indeed, from S-D equation at the 
Pisarski level approximation for low fermionic velocities, 
the self energy will acquire an imaginary part which can 
be interpreted as leading to a confinement 2 - for fermions 
of the theory. This instantly calls into question the va- 
lidity of using a pure plane wave type solutions for the 
computation of the scattering process. 

In this approximation, the nature of the solution 
changes considerably at the point where Va = 1- That is 
so because the radial integral g ets a log arithmic contribu- 
tion that is proportional to yl — v\ which overwhelms 
the leading contribution at the isotropic point. Alter- 
natively, in more physical terms, if photons and mass- 
less fermions move at precisely the same speed this is 

I 



By simple rescaling t' = r, x = Vpx 1 ', y — vpy' , a' T , — 
a T , a! x , — vpa x , and a' y , — Vpa y we can transform this 
theory into a new theory in which the fermionic part 

I 



1 (da' x , 
2e 2 v% \ dy' 



Thus, the effect of having vp ^ c s reduces to 
anisotropic couplings in the Maxwellian self-action of 
the gauge field. We denote these couplings as e 2 = 
e 2 Vp and e\ = e 2 v F . These two couplings can be in- 



"infinitely" different than having the photons that move 
faster than fermions. In the latter case it is natural to 
expect an "overscreening" behavior, in which constant 
exchange of fast photons ultimately leads to confinement. 
On the opposite side, with fermions moving faster than 
photons, we expect photons to be less effective in screen- 
ing the fermions, and thus less effective in generating 
their mass. We have also checked that in the isotropic 
limit the Pisarski's answer obtains and thus our results 
are not an artifact of the parametrization used. 

So far, we have shown that mass generation has a non- 
universal behavior, which arise due to the breaking of 
Lorentz invariance. Thus, it is natural for a cautious 
reader to wonder if, once N > N c , the renormalized ef- 
fective low energy theory is in fact Lorentz invariant or 
not. To begin with we emphasize that the fact that a 
flows to one by itself does not guarantee that the full 
Lorentz invariance will emerge unless both c s and vp are 
set equal; this in effect acts as an extra constraint. This 
is shown in Fig. ([7]) and in Fig. ([5]) where it is easy to see 
that even thought a — > 1 neither vp nor va converge to 
c s = 1 unless the above mentioned extra constaint is im- 
posed. To show that the full Lorentz invariance is indeed 
restored we must prove that both fermionic velocities vp 
and va flow to c s = 1 independently 

It is easy to see that Lorentz invariance can still be 
broken even with the fermionic anisotropy set to unity. 
Let us assume that the bare values of the Fermi and gap 
velocity are equal to each other but different from the 
gauge field velocity, i.e. vf = va, but vp ^ c s . The 
Lagrangian for this simplified theory i a 24 i 25 



(25) 

I 

of the action remains fully isotropic but an anisotropic 
Maxwellian term appears: 



I 

terpreted as two anisotropic charges. Such anisotropic 
charges can change the value of the critical number of 
fermions in the original theory, as already shown in Eq. 
(24). In contrast, within the isotropic QED3, the criti- 



C = ipi {jo{d T + ia T ) + r yiv F (d x + ia x ) + j2VF{d y + ia y )} ipi + (1 <-> 2) + 
1 J / da x da y \ 2 ( da T da x \ 2 / da T da y \ 2 1 
2e 2 1 \ dy dx J \dx dr J \ dy dr J J 




da' y , 
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cal number of fermions does not depend on e 2 . This is 
true as long as we have only one coupling constant, but 
once we introduce two different couplings this pleasing 
behavior is lost. This simple example shows that we can 
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FIG. 7: The RG-/3 function for vf (in arbitrary units) com- 
puted using the same method as in Ref. 7 but without setting 
da = c s = 1. This plot shows the parametric dependence of 
the fixed point on the value at which va is initially set in this 
calculation. 



restore Lorentz invariancc in the fermionic part of the 
action at the expense of breaking the Lorentz invariance 
of the Maxwellian part. 

To summarize, the above discussion shows that at a 
bare level there are two intrinsic anisotropies in this prob- 
lem; a is the fermionic one, while 9 = (e T /ej_) is the 
Maxwellian one. e T and e± are the couplings in the tem- 
poral and spatial directions respectively. Thus at a bare 
level of the QED3 theory relevant for the cuprates there 
are four coupling constants vf, va, e T , and ej_. We will 
now show that, even though they contain anisotropy both 
in the fermionic and Maxwellian terms, in the large N 
limit the relativity is ultimately restored, without any 
assumptions about the size of the anisotropy. 

To make good on the above claim notice that the ef- 
fect of fermions on photons is still described by Eq. (8). 
However, the gauge field stiffness is now: 



n (°) = — 

fJ,V 



(27) 



where the anisotropic Levi-Civita symbol is defined as 

^T7/jU — ^rr/fi, ^xr]fi — Q^xrifii ^Wl 1 — UVt^' ThuS, the 

effective Lagrangian of the theory can be written as: 

C A mso = \ (n{S + IV) a„{k)a„{-k) (28) 

where the effect of fermions has been introduced through 
the polarization function. This expression allow us to 
find the renormalizcd couplings, comparing the original 
bare gauge field stiffness with the screened one: 



r\(R) - n (°) + TT 



(29) 
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FIG. 8: The corresponding RG-flow of df and «a when charge 
renormalization is not taken into account. A line of fixed 
points at dp = va is apparent. 



We find that to the lowest order inl/N the renormalizcd 
couplings are: 



' 1 N 




A, 








4, 





N 
16k 
N 
16 



vfva 



(30) 



1 VF 



1 



r A + fc(2) Vf ) (31) 



where 1/k = 1/fcW + l/k^ 2 \ and fcW = 

The above one loop renormalization of the anisotropic 
charges allows us to set up the renormalization group 
(RG) equations for the beta-functions describing the flow 
of different couplings. These equations are rather com- 
plicated and we have been able to fully solve them only 
numerically. However, the following result is rather sim- 
ple and can be extracted in an analytic form: on general 
grounds we expect that the non trivial infra-red fixed 
point should remain once we introduce the anisotropy, 
even though its position in parameter space may change. 
In order to find the value of the renormalized couplings 
at the fixed point, we analyze the difference between the 
renormalized couplings of the Maxwellian action 













e -L/ R 







N 
16 



Vp 



+ 
1 



(32) 



v A v F - 



VF 



If we now rearrange our RG equations so as to focus on 
the beta-function for this difference between the renor- 
malized charges, f3 eT _ e± , the above equation implies that 



N 
16 



fed) 



fc(2) 



+ 



(...) 



,(33) 
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where (• • •) can be rewritten in terms of beta- functions for 
all other couplings and thus must vanish at the putative 
fixed point. Clearly, noting that N is an arbitrarily large 
number, and that k^ 1 ' share the same sign, it follows that: 



(R) 

JR) (R) _ Vp_ 
v F v A , R) 

V A 



= 



(34) 



from where — 1, and, given that the theory is fully 
invariant under the exchange vf <-» Va, it follows that 

( R) 

v F ' = 1. Putting this information back in the flow equa- 
tions it is clear that 1/e^ and 1/e^ themselves diverge 
with the same slope at the fixed point, and therefore their 
ratio 9 — (e T /e±) — > 1. This result shows that Lorentz 
invariance is restored and thus the previous results^i 
remain valid. However, we have made it clear that the 
physics behind the restoration of full Lorentz invariance 
follows a path more subtle than previously explored: at 
the infra-red fixed point the relativity is restored due 
to the interplay between the velocity and charge renor- 
malizations, the velocity renormalization by itself being 
insufficient to fully restore relativity of the theory. 



III. CONCLUSIONS 

We have proposed a simple criterion that allows an ex- 
plicit computation of the critical number of fermions N c 
in a theory that contains intrinsic anisotropics. We have 
checked, at the Pisarski's level approximation, that this 
criterion captures the functional dependence of N c (va) ~ 
1/va in the case in which an explicit expression can be 
obtained from analytic calculations. 

Our criterion suggests that lattice simulations should 
be performed in such a way that important symmetries of 
the theory, namely vp <-> va, are protected. Otherwise, 



there is a danger of obtaining spurious results. In lattice 
QED 3 , it seems worthwhile to investigate the existence 
of a possible confined phase in the region of the param- 
eter space where fermionic velocities are small compared 
with the gauge field velocity. Another venue that remains 
to be explored, is the possible usefulness of similar cri- 
teria for the analysis of CSB or other non perturbative 
phenomena in other physical systems that also feature 
anisotropic couplings. 

Finally, we have shown that the velocity anisotropy in 
(2 + 1)QED does affect the number of critical fermion 
flavors N c at which chiral symmetry is broken due to the 
phenomenon of mass generation, even though the large 
N theory remains fully relativistic in its critical phase. 
Surprisingly, N c is a non-monotonic function of vp and 
i>A, and, depending on the specific value of the ratio 
vf/c s , different regimes emerge. Our results imply that 
if phase fluctuations destroy the superconducting order 
in underdoped cuprates, we should expect a protected 
chirally symmetric critical phase - i.e. the pseudogap 
within this theory - as doping decreases, before we reach 
the antiferromagnetic region in the phase diagram. De- 
tails of how will this sequence take place depend on the 
specific value of the gauge field velocity c s for different 
compounds. We hope that our results will contribute 
to better understanding of the quantitative issues that 
surround the value of -/V c in various effective theories and 
motivate further research on the anisotropic incarnations 
of the QED 3 theory. 
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In Ref.[24] the authors claim that fermionic anisotropy 
goes to one, and thus they obtain the simplified model we 
used to show that this can be thought as a model in which 
there are anisotropic couplings in the Maxwellian term. 



Given that the number of fermions in physical systems is 
modest it would be interesting to analyze the consequences 
of having such anisotropic couplings. 



